Colloidal glass transition observed in confinement 
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We study a colloidal suspension confined between two quasi-parallel walls as a model system 
for glass transitions in confined geometries. The suspension is a mixture of two particle sizes to 
prevent wall-induced crystallization. We use confocal microscopy to directly observe the motion of 
colloidal particles. This motion is slower in confinement, thus producing glassy behavior in a sample 
which is a liquid in an unconfined geometry. For higher volume fraction samples (closer to the glass 
transition), the onset of confinement effects occurs at larger length scales. 
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Glasses are typically formed by rapidly quenching the 
temperature of a liquid, resulting in an amorphous liquid- 
like microstructure with macroscopic solid-like behavior. 
Upon approaching the glass transition, the temperature 
might be changed by only a factor of two while simultane- 
ously the viscosity of the liquid grows by many orders of 
magnitude [l[ . A conceptual microscopic explanation for 
the viscosity growth is the idea of dynamic length scales: 
in order for molecules in the material to rearrange, they 
must move together as a group. As the glass transition 
is approached, the increasing size of these groups relates 
to the increasing macroscopic viscosity P, 0, H,H, d, [f| . 

An important way to probe these length scales is to 
study the behavior of glass-forming systems when they 
are confined, to constrict the range of accessible length 
scales @, 0, H, I, M El [H, El M- Intriguingly, the 
glass transition temperature T a can be both smaller or 
larger in confined geometries [ill, E2, EHj even for the 
same material 0, 0, fl3 |. Experiments and simulations 
suggest that the interaction between the confining sur- 
face and the sample is crucial. For strong interactions (or 
atomically rough surfaces) the glass transition happens 
"sooner," that is, confinement increases T g by slowing 
motion near the surfaces 0, 0, H, EH E3 ■ Likewise, for 
systems that weakly interact with the walls, T g is typi- 



cally smaller 0, 0, [Tl|. However, a clear explanation of 
these phenomena is still lacking. As it is difficult to get 
details out of experiments [2j , the use of computer simu- 
lations to visualize the motion is important @, [1, [H, E3] ■ 

We use confocal microscopy to directly visualize the 
motion of colloidal particles, which serve as a model sys- 
tem for the glass transition in confinement. Colloids un- 
dergo a glass transition in bulk samples as the solid par- 
ticle volume fraction <f> is increased |, SEE [3. At 
high volume fraction near the colloidal glass transition 
{4> g « 0.58), particles move in rearranging groups charac- 
terized by a length scale of ~3-6 particle diameters [1, Hj], 
similar to simulations [101 ] . In this manuscript we study 
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a mixture of two sizes of colloidal particles confined be- 
tween two quasi-parallel plates, with a plate gap as small 
as 3.0 large-particle diameters. In our system confine- 
ment induces the glass transition "sooner," at concen- 
trations for which the bulk behavior is still liquid-like. 
Studying the glass transition in confinement may help us 
understand the glass transition in the bulk 0. Further- 
more, understanding the properties of confined liquids 
has relevance for lubrication [3, E3 > dusty plasmas [l|[ , 
and the flow of glassy complex fluids through microfluidic 
devices [l7j . 

Our colloidal samples are poly-methyl(methacrylate) 
articles, sterically stabilized to prevent aggregation 
"l5|. We use a mixture of two particle sizes, with radii 
flsmaii = 1-18 /im and ai argc = 1.55 /mi. While the parti- 
cle polydispersity is low (~5%), the mean particle radii 
are only known to within ±0.02 /mi. The mixture of two 
particle sizes prevents crystallization which would other- 
wise be induced by the walls [1, EH EH- In each sample, 
the small particles arc dyed with rhodamine dye, and the 
large particles are undyed. We use a mixture of cyclo- 
hexylbromide and decalin as our solvent, to match the 
density and index of refraction of the particles; the parti- 
cles are slightly charged in this solvent 20] . The viscosity 
of the solvent is 2.25 mPa-s. We examine four different 
samples A-D, with properties listed in Table I. 

We observe our samples using confocal microscopy (20l . 
2lj . As the larger particles are not dyed, we only see 
the smaller particles. We use a fast confocal microscope 
(VT-Eye from Visitech, International) with a 63 x air 
objective (N.A. 0.70) to scan a volume 50 x 50 x 20 /im 3 
once every 2.0 s over a period of an hour. We analyze the 
images offline to locate the positions of visible (smaller) 
particles, with a resolution of 0.05 /mi in x and y (parallel 
to the walls) and a resolution of 0.1 /jm in z (parallel to 
the optical axis) . We then track their motion in 3D [2(| • 
Our sample chambers are made by placing a micro- 
scope coverslip at a slight angle, supported by a small 
piece of mylar film (thickness 100 /im) at one end and 
resting directly on the microscope slide at the other end 
[T9L |22| . The ends and sides are sealed with UV-curing 
epoxy. Thus a thin wedge-shaped chamber is formed with 
an angle ~ 0.4°, ensuring that locally the walls are essen- 
tially parallel, and allowing us to study a single sample 
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TABLE I: Characteristics of the four samples studied. The 
number ratio -/V ama n/iVi ar g is determined by counting parti- 
cles in several fields of view using DIC (differential interfer- 
ence contrast) microscopy. The total volume fraction (j> tot is 
determined using confocal microscopy, by counting the num- 
ber of small particles seen in a given imaging volume, using 
the known number ratio to determine the number of large 
particles present, and then using the particle sizes and the 
imaging volume size to compute 4>tot- Additionally (f> to t and 
A'smaii/A'i ar g e was confirmed in samples B-D by direct 3D con- 
focal microscopy observation, where the particle sizes could 
be easily distinguished and counted; the results were in agree- 
ment with the DIC measurements. The volume fractions 
of the small species and large species, <f> s and 4>i, are cal- 
culated from the other two quantities. The uncertainties of 
A'smaii/A'iarge are ±5%, and the uncertainties of (f>tot are ±8%. 
In particular, note that samples A and C likely do not have 
the same volume fraction, but it is unclear which has the 
larger <f)tot- Samples B, C, and D are prepared by dilutions of 
one stock sample and thus all have the same A'smaii/A'iarge- 



Sample 


N small / -^largc 


4>s 


4>i 


4>tot 


A 


3.5 


0.26 


0.16 


0.42 


B 


3.0 


0.13 


0.10 


0.23 


C 


3.0 


0.24 


0.18 


0.42 


D 


3.0 


0.26 


0.20 


0.46 



at a variety of different confinement thicknesses [13, |22| • 
The glass surfaces of the coverslip and slide are un- 
treated. In experiments with sample A, we find that some 
colloidal particles stick to these surfaces. The stuck par- 
ticle coverage is typically 10% - 20% of the area. In a 
second series of experiments done with samples B-D, no 
particles were stuck. Reassuringly, we find little depen- 
dence of the behavior on the number of stuck particles in 
the results discussed below [231 ] . 

For sample A, measuring the positions of the stuck par- 
ticles allows us to accurately measure the sample thick- 
ness. While the uncertainty in locating individual par- 
ticle positions in z is 0.1 /im, by averaging data from 
tens of stuck particles over hundreds of images we locate 
their mean z position to better than 0.005 /im. Thus 
the effective thickness H of each experimental data set 
is determined to within 0.01 /im, and is the range in z 
available to the centers of the visible particles. In this 
manuscript our thicknesses are reported in terms of H. 
The true surface-to-surface thickness of a sample cham- 
ber is found by adding 2et sma n = 2.36 fim to H. 

For the first series of experiments, we study the be- 
havior of sample A (<p ks 0.42) as a function of thick- 
ness. We quantify the particle motion by calculating 
the mean square displacement (MSD), (Air 2 ) = ((xi(t + 
At) - x,(t)) 2 ), where the average is taken over all par- 
ticles i and all initial times t, and a similar formula ap- 
plies for {Ay 2 ) and (Az 2 ). We find that (Ax 2 ) « (Ay 2 ) 
for all our experiments; we report our results for the x 
direction, the direction over which the sample chamber 
has constant thickness. We first consider the results for 
motion parallel to the confining plates, (Ax 2 ), shown in 
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FIG. 1: Mean square displacements, (a) Data for sample 
A, showing motion parallel to the walls, for thicknesses H = 
bulk, 16.28 /Km, 11.06 ^m, 9.41 /jm, and 6.92 /im (from top 
to bottom). The dashed line has a slope of 1 and indicates 
the expected motion for a very dilute bulk suspension of par- 
ticles with radius a sma ii. (b) Similar to (a), but for motion 
perpendicular to the walls. Data are ordered by thickness as 
At — > co, as in (a). 

Fig. [IK a) . The upper bold line shows motion in an un- 
confincd region and is reproducible for all chamber thick- 
nesses H > 20 /xm. For this sample, the motion in the 
unconfined region is nearly diffusive, with the MSD grow- 
ing almost with slope 1 on the log- log plot. This behavior 
is similar to monodisperse samples with a volume frac- 
tion of </> < 0.4 4]. In other words, this sample is far 
from the glass transition, <fi g ~ 0.6 [m [24| . In thinner 
regions, the motion slows, as seen in the sequence of solid 
curves below the top bold curve in Fig.QJa). This slowing 
starts at a thickness of H ~ 16 /urn (2nd curve from top) 
and slows dramatically for thinner samples; note that 
Fig. QJa) shows a log-log plot and thus for the thinnest 
region shown (bottom curve, H = 6.92 /im), to move 
a distance (x 2 ) = (a sma n/3) 2 it takes a time scale 180 
times larger than for the bulk region data (At = 500 s as 
compared to 2.8 s). 

These results suggest that confinement induces glassy 
behavior, with the influence of confinement beginning at 
H w 16 /im rs 14a sma u w 10ai arge for this sample. For 
the lower curves in Fig. [lja), the characteristic behavior 
of a "super-cooled" sample is seen: at shorter lag times 
(At < 100 s), the MSD has a plateau, while at longer 
lag times, the MSD begins to rise again 0,0]- (For short 
time scales, particles are trapped in cages formed by their 
neighbors, causing the plateau in (x 2 ). At longer time 
scales, these cages rearrange 0, For the thinnest 

region (bottom curve), the particles remain localized for 
the duration of the experiment. 

The slowing is also seen in motion perpendicular to the 
walls, quantified by (Az 2 ), shown in Fig.QJb). Moreover, 
in comparison with Fig. [IJa), it is seen that the motion 
perpendicular to the walls, (Az 2 ), is slowed even more 
so than motion parallel to the walls, (Ax 2 ). This is not 
surprising, given that particles close to the walls cannot 
move toward the walls at all, whereas motion parallel to 
the walls is less restricted. 

More than merely constricting motion, the walls also 
induce a layering of particles, as seen in Fig.(2ja), similar 
to simulations [8|, [25j . The layering is most pronounced 
immediately adjacent to the walls. The centers of these 
peaks are not at the precise distance a sma ii from the walls, 
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FIG. 2: (a) Particle number density n sma ii(.z) as a func- 
tion of distance z across the sample cell. Additional parti- 
cles are stuck to the walls of the sample cell (not shown in 
the plot) which have centers located at z = 0.00 /im and 
z — H — 11.06 /im. This data corresponds to the middle 
curve in Fig. [Ha), that is, a sample with moderately slowed 
dynamics, (b) Mean square displacement parallel to the walls 
(x) and perpendicular to the walls (z), as a function of z. The 
displacements are calculated using At = 100 s. The vertical 
dotted lines indicate the positions of the peaks from part (a). 
Data shown are from sample A ((/> « 0.42). 



but are slightly offset toward the interior of the sample. 
(The centers of the stuck particles indicate the maximum 
possible extent in z that particles could be located, and 
correspond to the "feet" of the data shown in Fig. [2ja) 
at z = and 11.06 /^m. These particles are not counted 
in n sman shown in Fig.[^a).) 

The layers influence the dynamics, as seen in Fig.[2jb), 
which shows how (Ax 2 ) and (Az 2 ) depend on z. The 
displacements are calculated using At = 100 s, as a rep- 
resentative time scale over which particles begin to move 
out of this cage, although the results do not depend on 
this choice and are similar for caged behavior (At = 2.0 s 
for example). Particles in the layers [the peaks of n(z)] 
have smaller vertical displacements, as seen by the dips 
in (Az 2 ) (heavy line). The implication is that particles 
in layers are in a preferred structure and less likely to 
move elsewhere [8Ll25j|. 

Surprisingly, the layers do not appear to influence the 
motion parallel to the walls, as seen by (Ax 2 ) (thin line), 
which does not depend on z. (The slight dip in (Ax 2 ) 
seen at the largest values of z is not reproducible in other 
data sets.) This seems counterintuitive as hydrodynamic 
interactions with the wall normally result in reduced mo- 
tion for nearby particles [26| . We speculate that the cage 
dynamics dominate particle motion, rather than hydro- 
dynamic influences [27] . For example, if a particle is 
pulled by an external force in a direction parallel to the 
walls, other particles would be forced to rearrange, which 
is probably the most significant contribution to the drag. 



Particle rearrangements would be even more constrained 
for a particle pulled perpendicular to the wall, thus ex- 
plaining why we observe slower z motion [27| • Simply 
put, the high volume fraction likely results in hydrody- 
namic screening. 

Thus while confinement causes the layering of particles 
near the walls, this layering does not appear directly re- 
sponsible for the slowing of the particle motion. Rather, 
the layering seems to be an additional influence on the 
motion in the direction perpendicular to the walls, as 
seen in Fig. [Hb), but only a minor influence compared 
to the overall fact of confinement. Note that results do 
not appear to depend on having an integral number of 
well-defined layers between the walls Q . The overall dy- 
namics slow smoothly and monotonically as the confining 
dimension decreases. 

Our observation that the layers closest to the wall have 
slower motion perpendicular to the walls agrees qualita- 
tively with previous experiments (l2l . [l3l [l4| and sim- 
ulations 0] which suggested that surface layers may be 
glassier than the interior. However, we note in our ex- 
periment this is strongly directionally dependent. The 
slowing is most easily seen if (Az 2 ) can be measured in- 
dependently of the other two directions. 

As noted earlier, the growth of dynamic length scales 
has been observed as the glass transition is approached 
in a bulk material [H, 0, S 0, H @|- F° r our colloidal 
samples, this implies that samples with a larger <j) should 
exhibit stronger confinement effects. To check this, we 
took data from samples B, C, and D at various thickness. 
Qualitatively the data resemble that shown in Fig. [lja). 
To capture the H dependence, Fig. [3] shows the values 
of (Aa; 2 ), at fixed At = 100 s, as a function of H for 
the different samples. Consider the solid triangles, corre- 
sponding to sample D. For H > 50 [im, (Ax 2 ) is essen- 
tially constant. At H < 50 /im, the data start showing 
a strong H dependence, suggesting a confinement length 
scale of H* « 50 /xm. For the solid symbols, an increase 
in H* is seen as <f> increases, from approximately 10 /xm 
to 50 /im, confirming that there is a growing length scale 
as the glass transition is approached. These length scales 
are significantly larger than those seen for dynamical het- 
erogeneities in monodisperse samples, which arc 4 8 /mi 
Q. However, this agrees with simulations which found 
a confinement length scale significantly larger than the 
mobile cluster size Q. In Fig.[3j sample A has a smaller 
value of H* relative to sample C, which may be due to 
the excess of small particles in sample A; see Table I. 

We find that confinement slows the motion of colloidal 
particles and thus induces a glass transition to occur 
sooner than normal, in other words, at volume fractions 
for which the bulk behavior is liquid-like. Simulations 
suggest the roughness of the walls is crucial to this slow- 
ing and we plan to vary this in future experiments. 
However, we note that our data show slowing both with 
completely smooth walls (samples B, C, and D) and walls 
with isolated stuck particles (sample A). In contrast to 
our work, rough walls in simulations are composed of 
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FIG. 3: Value of (Ax 2 ) at At = 100 s, as a function of thick- 
ness H , for samples with <j> as indicated. The open circles cor- 
respond to sample A with AT slna u/iVi ar g e = 3.5, while the solid 
symbols correspond to samples B-D with iV sma ii /Marge = 3.0. 
The lines are drawn to guide the eye. The plateau for each 
data set indicates behavior corresponding to the bulk, whereas 
the downturn at low H gives an idea of the length scale at 
which confinement becomes important. 



particles fixed in a liquid- like structure @, [f| ■ This pre- 
vents layering of adjacent particles and restricts motion 
parallel to the walls. Thus the glass transition in con- 
fined samples occurs sooner (at higher temperatures Q 
or lower densities |8|). In our experiments, particle mo- 
tion parallel to the wall is not noticeably inhibited, as 
seen in Fig. [5J Yet, we still find the glassy behavior oc- 
curs sooner: at constant volume fraction, the dynamics 
are slower as the confining dimension decreases. Thus it 
seems that the important effect in our experiments is sim- 
ply the restriction of motion perpendicular to the wall, 
close to the surface of the wall. 
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